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We analyze the nonlinear active microrheology of dense colloidal suspensions using a schematic model of mode-coupling theory. 
The model describes the strongly nonlinear behavior of the microscopic friction coefficient as a function of applied external force 
in terms of a derealization transition. To probe this regime, we have performed Brownian dynamics simulations of a system of 
quasi-hard spheres. We also analyze experimental data on hard- sphere-like colloidal suspensions [Habdas et al, Europhys. Lett., 
2004, 67, 477]. The behavior at very large forces is addressed specifically. 



1 Introduction 

Microrheology is a developing technique to acquire local in- 
formation on the viscous and elastic properties of complex 
fluids and soft matter.'^ In active microrheology one uses 
a mesoscopic colloidal probe particle manipulated by an ex- 
ternal driving mechanism such as laser tweezers or magnetic 
forces, and measures the response of that particle to the driv- 
ing. The two major complications in making sense of the ob- 
tained data are, firstly, that one usually deals with suspending 
host liquids that have mesoscopic structure on length scales 
comparable to the size of the probe, calling into question sim- 
plifying assumptions treating the host liquid as a continuum. 5 
j Second, even moderate external driving is sufficient to enter 
the nonlinear-response regime of soft host liquids. On the 
other hand, if one is able to understand the implications of this 
situation, the technique in principle gives access to detailed in- 
formation about structure-dynamics relationships of complex 
soft matter. This makes active microrheology an ideal tool to 
probe, e.g., cellular environments .DEEl 

Here, we focus on the application of constant-forcing ac- 
tive microrheology to colloidal glass formers, different from 
variable forcings such as parabolic traps moving at cons tant 
velocity, 8 or where the tracer moves at a constant speed.^ESl 
Approaching the glass transition, the dynamics of the host liq- 
uid becomes increasingly slow, heterogeneous, and stretched 
over many orders of magnitude in time. Microrheology might, 
if properly understood, provide valuable insights into the mi- 
croscopic origins of this slow structural relaxation dynam- 

a Fachbereich Physik, Universitdt Konstanz, 78457 Konstanz, Germany. 

b Departamento de Fisica Aplicada, Universidad de Almeria, 04.120 Almeria, 

Spain. E-mail: apuertas@ual.es 

Institut fiir Materialphysik im Weltraum, Deutsches Zentrum fur Luft- und 
Raumfahrt (DLR), 51170 Koln, Germany. 

d Zukunftskolleg, Universitdt Konstanz, 78457 Konstanz, Germany. 

t Present address: Leibnitz Institute of Polymer Research Dresden, Hohe 

Strafie 6, 01069 Dresden, Germany 



ics.EE^The slow relaxation processes render highly nonlin- 
ear the relation between the measured steady- state velocity of 
the probe and the externally applied force. A pronounced su- 
perlinear rise in the velocity is observed for forces far exceed- 
ing the scale set by thermal energy and the particle size This 
force threshold can be interpreted as the strength of nearest- 
neighbor cages that are broken by strong forcing before they 
can relax due to structural relaxation. 

We have recently presented a microscopic theory of nonlin- 
ear force-driven active microrheology for dense colloidal sus- 
pensions:^ based on an integration-through transients (ITT) 
framework and mode-coupling approximations, equations for 
the microscopic friction coefficient ^(F QX ) were obtained, de- 
fined by the steady-state relationship 



C(v), 



(1) 



where angular brackets denote the steady- state ensemble av- 
erage reached at long times. These equations have a fully mi- 
croscopic foundation as they are based on the Smoluchowski 
equation for the colloidal system (neglecting hydrodynamic 
interactions). They only require the equilibrium static struc- 
ture factor as input. This, however, makes them also rather 
complicated to solve. In order to understand generic features 
of the equations, we have devised schematic models, i.e., ad- 
hoc simplifications of the original equations that are much eas- 
ier to solve. The simplification essentially amounts to drop- 
ping spatial information (wave- vector dependences), at the 
cost of introducing a limited number of fitting parameters that 
replace the static- structure-factor input. One such schematic 
model was already presented earlier,^ but this model did not 
capture some features of the high-force limit of the dynam- 
ics observed in experiment and simulation. In the following, 
we present an extended model that levies this limitation by 
taking into account separately the symmetry-breaking direc- 
tion along the external force, and directions perpendicular to 
this. We are thus able to discuss the qualitatively different be- 
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havior of probe fluctuations along those directions, and their 
effects on the friction in the steady- state probe velocity. The 
model is checked with computer simulations of a Brownian 
system, which we perform, and experiments taken from Hab- 
das et alM^The model correctly describes the qualitative fea- 
tures of the friction coefficient as a function of density and 
external forcing, and reproduces them quantitatively with a 
sensible choice of the parameters for both the simulations and 
experiments. The tracer position correlation function, a key 
quantity in the theory, behaves similarly in simulation and the- 
ory. 

2 Computer Simulation 

We performed molecular-dynamics simulations for quasi- 
hard-sphere particles governed by a Langevin equation, in the 
following referred to as Brownian dynamics (BD) simulations 
as hydrodynamic interactions are neglected. The system is 
poly disperse to avoid crystallization (flat distribution of radii, 
with half- width 10% of the average radius), and the interac- 
tion potential is given by V(r) = ksT (r / a) -36 , where a is the 
particles' center- to-center distance. ksT = 1 and the average 
radius of a host-liquid particle, a = 1 , set the units of energy 
and length. All particles have the same mass, m—\. The sol- 
vent friction coefficient, which also sets the Langevin forces, 
is fixed to = 50 in these units. The near-equilibrium dynam- 
ics of this system has been analyzed in detail beforej 15 ! 16 ! es- 
tablishing a mode-coupling glass transition at an overall pack- 
ing fraction (p « 0.595. At the highest packing fraction in- 
vestigated here, (p = 0.62 the system cannot be equilibrated 
in the available computer time; here we report measurements 
obtained after a waiting time t w = 2.5 x 10 4 after setting up 
the simulation runs. For forces exceeding F « 35ksT/a, this 
is long enough to obtain results that do not depend on t w \ for 
smaller forces, aging effects may still affect the data at this 
highest packing fraction. 

The simulation box is 8 times longer in the direction of the 
applied force than in the transversal directions, to allow the 
study of longer trajectories of the tracer, and contains 1000 
particles, as shown in the snapshot in Fig. [T] At the highest 
density studied here, (p = 0.62, the width of the box is « 9.5a, 
enough to relax the structural effects caused by the moving 
probe. At time t = 0, one particle is selected at random to be 
the probe particle (hence, the probe radius a s = a = 1), and 
a constant external force F ex = (F ex ,0,0) is exerted over it 
(in addition to the interparticle, friction and Brownian forces), 
and its trajectory is monitored. When the tracer has travelled 
a distance equal to half the box length, a new probe is selected 
from scratch, and the average velocity is obtained from the av- 
erage trajectory. The simulation scheme is identical to the one 
used recently. 14 For small forces, trajectories were recorded 
up to t = 10000, before selecting a new probe. In these cases, 




Fig. 1 Snapshot of the simulated system at packing fraction 
(p = 0.55, showing the box geometry and an examplary run with a 
randomly selected probe. An arrow indicates the direction of the 
applied force. 

sampling was performed over 300 trajectories obtained from 
different initial conditions and tracer particles. 

From the simulation runs we thus obtain the probe re- 
sponse in terms of the friction coefficient ^(F QX ) at various 
packing fractions, using the steady-state relationship, Eq. ([]]). 
We also analyze the dynamics of density fluctuations, char- 
acterized by the tagged-particle density correlation functions 
<I>|W = (P s *(^t)p s (q)). Here, p s (q) = exp[iqr s ] is the tagged- 
particle density fluctuation to wave vector q. In an isotropic 
homogeneous liquid, the correlation function depends on the 
wave vector only through q = \q\; however, in the present 
context, application of the external force degrades this spher- 
ical symmetry to a merely rotational one. As above, angular 
brackets denote ensemble averaging; here, however we distin- 
guish equilibrium from steady-state averages, leading to, re- 
spectively, transient and stationary correlation functions. 

3 Schematic Models 

In this section, we summarize the main equations defining 
the schematic models. Mode-coupling theory (MCT) in the 
integration-through transients framework uses temporally and 
spatially resolved transient density correlation functions to de- 
scribe the nonlinear response of the system to the strong ex- 
ternal perturbation. In the case of constant-forcing active mi- 
crorheology, the transient tagged-particle density correlation 
function is the basic dynamical quantity. Acknowledg- 
ing the symmetry-breaking direction of the external force, it 
will be useful to distinguish for the probe particle the direc- 
tions q || F ex and q JL F ex . If q has a component parallel to 
the force, the displacement of the particle due to the external 
force causes the correlator to become complex valued, as it is 
the Fourier transform of a tagged-particle density profile that 
is not centered on the origin. Formally, this is a consequence 
of the non-Hermitian time-evolution operator in the Smolu- 
chowski equation without detailed balance. If q is perpendic- 
ular to F ex , the correlator stays real. The internal forces acting 
on the tracer are determined by tracer- and host-liquid den- 
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sity fluctuations. Whitin MCT, the latter are captured by the 
intermediate scattering function, § q (t), which, in an appropri- 
ate thermodynamic limit can be taken as the equilibrium one 
(being real and exhibiting full rotational symmetry). 

The schematic model describes the dynamics of the system 
through wave-number independent functions §(t) and § s a (t) 
that serve as proxies for the full correlation functions; here, 
a G { || , -L} tells apart a complex- and a real- valued tracer cor- 
relator. The neglect of spatial resolution is motivated by equi- 
librium MCT, where a factorization theorem states that all 
length scales are strongly coupled and arrest simultaneously 
at the MCT transition.^ The two kinds of tracer correlators 
capture the dominant effects of the external force: some cor- 
relators turn complex, while others remain (almost) real. The 
non-Markovian equations of motion analogous to the micro- 
scopic Mori-Zwanzig equations read 

3 / <|)(0+r|<|)(0 + ^m(r-O9/'<l>(O* / | =°> ( 2 ) 

a^ a (0+<{ra(0+^<(^-O«(O^} =o. O) 

Here, T models the short-time dynamics of the host-liquid par- 
ticles; we choose time units such that T = 1. The short- time 
dynamics of § s a (t) is affected by the non-Hermitian nature 
of the underlying time-evolution operator and is modeled by 
<a| = - /K F ex ) and co^ = F s . The coefficient F s could be 
used to describe a short-time diffusion coefficient of the probe 
that differs from that of the host-liquid particles; we set T s = 1 
for simplicity. The strength of the external force relative to 
internal ones at short times is measured by Ko. 

The memory kernels m(t) and m s a (t) are approximated in 
MCT as nonlinear functionals of the density correlators them- 
selves. We take for the schematic model 



m(t)=v l §(t)+v 2 §(t) 2 



(4) 



m\it) = [v s ^f (t)m+W ± (tMt)\ /(l-^ ex ), (5) 

<(0= [v!^(0*(^) +^*f|(^)<l>(0] / C 1 + (kj.^™) 2 ) • 

(6) 

All coupling parameters v, and v • are taken to be real and posi- 
tive; they model force-free equilibrium structural correlations. 
Equation ^ just specifies the well-known F12 model that is 
regularly used to analyze linear-response dynamics of glass 
forming liquids within mode-coupling theory.^ The quadratic 
polynomial mimics the feedback mechanism termed 'cage ef- 
fect' that causes the slowing down in the structural relaxation 
of the host liquid. The bilinear coupling to § s a (t) and §(t) 
in Eqs. ([5]) and ([6]) mimics that internal forces on the probe 
relax via its motion and via rearrangements of the surround- 
ing particles. Both tracer motion parallel and perpendicu- 
lar to the external force can relax the local friction. Spatial 



inversion symmetry demands that (|>jj(f) and §**(t) are cou- 
pled, which is assured by the complex conjugate appearing 
in mjj(f). The ensuing force driven derealization transition 
can be mapped out using the external force as single relevant 
control parameter.^! Since m s ± (t) is real- valued, only the real 
part 9t(|)jj (t) can enter in Eq. ([6]). The external force enters 
differently in mjj(f) and m s ± (t): the kernel mjj(f) is complex- 
valued, and F ex suppresses both parts, while m s ± (t) needs to 
stay real-valued. Including F ex in the denominator ensures 
the correct physical behavior that large forces decrease the 
coupling strength. Equations §5§ and ^ explicitly express 
the symmetry under inversion of the force: (|)| (— F) — (|>jj*(F) 
but § s ± (— F) — ^(i 7 ). The positive parameters Ky and 
measure the effective force in the directions dominantly paral- 
lel and dominantly perpendicular to the external force. They 
should, by analogy to the microscopic model, be smooth func- 
tions of the thermodynamic control variables,^ and we take 
for simplicity Ky = Kq which was the choice in the original 
schematic model. 14 

Finally, the friction coefficient is expressed through dy- 
namical correlation functions via a nonequilibrium general- 
ization of the Green- Kubo relation central to the ITT frame- 
work.^ Following a mode-coupling approximation, we arrive 
at, schematically, £ = ^0 + A£, 

r-oo poo 

AC/Co = Ai||r J y A9tyjj(*)<K0 +^ r */ *<B_(0<K0, ( 7 ) 

where the parameters id a >0 replace angular-dependent cou- 
pling coefficients that are given by the equilibrium structure 
functions in the microscopic theory. Again this MCT approxi- 
mation expresses that the friction on the probe arises from the 
cumulated transient fluctuations of probe and host-fluid den- 
sities. From the correlation function a =||, only the real part 
enters due to spatial inversion symmetry. Equations ^ to ([7]) 
specify our schematic model completely. They can further be 
motivated by an ad-hoc simplification of the full microscopic 
MCT, restricting wave- vector integrals to four wave vectors, 
q || F ex and q _L F ex M^For each q also — q has to be taken to 
obey the required symmetry § s _^{t) = § s £(t)', this justifies the 
appearance of the complex conjugate in ([5]). 

The earlier schematic model for active microrheologyE2l 
only included a single complex correlator, concentrating on 
the novel transition that a particle is pulled mobile by a finite 
external force. The present three-correlator model reduces to 
this two-correlator model upon setting v| = ju± = 0. Also, 
for F QX = 0, we obtain <|>jj(f) = <|>j_(0 and the model reduces 
to the Sjogren model of tagged-particle dynamics in the near- 
equilibrium glass forming liquid, with a single probe-coupling 



parameter v s 



; the latter is a standard model for anal- 



ysis of linear-response dynamics of glass forming liquids, e.g., 
in terms of scattering spectra.^ 
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The present model additionally captures a nontrivial limit 
of the friction at high force, A^oo = A^(F QX — >> oo). The model 
with ju± = predicts = as the presence of the term oc 
iF ex in co| leads to arbitrarily fast oscillations in (|)| (t); in the 
integral determining A^ these lead to cancellations such that 
the overall integral vanishes. Only when ju± ^ does one get 
an additional contribution to A^ in the large force limit. If 
we furthermore neglect all memory kernels, we obtain as the 
low-density (weak-coupling) limit of the schematic model the 
ratio A^oo/A^i 7 ^ — >> 0) = /li±/(jli\\ +A / ±). This ratio is known 
to equal 1/2 exactly in the low-density limit^] checked also 
in our simulations. We thus, for simplicity, set m = ju± = /u 
throughout in the following. 

4 Data Analysis 

The schematic model presented above has several fit parame- 
ters. Two of them, v\ and V2, do not depend on the probe at 
all; they specify the dynamics of the host liquid and its vicin- 
ity to the glass transition which occurs for some critical cou- 
pling (vpV^). Thus, these two parameters can in principle be 
determined independently. In the cases we analyze, the only 
physical control parameter is the packing fraction of the sus- 
pension, (p. One demands that the fit parameters follow a lin- 
ear relationship (vi,V2) = (vpV^Xl+s) wither (cp — (p c )/cp c , 
where (p c is the glass-transition packing fraction of the suspen- 
sion. The critical point is then fixed by demanding the MCT 
exponents of the schematic model to match those found in ex- 
periment or simulation. In practice, this fixes the pair (vf, v^), 
and one fits £ for each packing fraction, biasing it to follow the 
stated linear relationship. 

Two parameters v\ and v s 2 describe the coupling strength of 
the probe to the host suspension; their sum v s could be deter- 
mined from analyzing linear-response measurements indepen- 
dently. In principle, they depend on the density of the host; we 
will in the following thus allow v s to increase with increasing 
packing fraction. Inspection of the microscopic model and its 
symmetries 18 suggests the fixed ratio v\/v s 2 = 2. The remain- 
ing parameters and Ky are global parameters (fixed differ- 
ently for the simulation and for the experiment to reflect the 
difference in systems) specifying the relative influence of the 
external force on the memory kernel. 

In our fits of both BD simulation and experimental data, we 
aimed at keeping as many parameters fixed as is physically 
plausible. First, we set v c 2 — 2, resulting in v\ — 2{\fl— 1), 
which renders the near-equilibrium asymptotic behavior of the 
schematic model close to that observed for the hard-sphere 
glass transition. Furthermore, e((p) is estimated from the 
known linear-response regime for the simulation data; for the 
fits to the experimental data we require that the obtained val- 
ues of 8 are of the same magnitude for (p not too close to the 





Fig. 2 Microscopic friction coefficients £ as a function of external 
force F ex : symbols are (a) Brownian dynamics (filled) and 
hard-sphere simulations of Ref. [20](open); (b) experiment of Ref. 
1T21 evaluted from the original velocity data (Fig. [3]) and translated to 
dimensionless quantities by setting = 0.095 mg/s and 
lpN = 271.9£#r la. Lines are fits with schematic models of 
mode-coupling theory, taking only fluctuations in direction of the 
force into account (dashed), or also in the perpendicular direction 
(solid lines). Horizontal dot-dashed bars in (a) give the quiescent 
self-diffusion coefficients. 



critical point. To capture the density dependence of the high- 
force plateau, we also increase /n with increasing density. 

We first turn to a discussion of the friction coefficient 
shown in Fig. [2] for our Brownian dynamics simulations, and 
for the experiment by Habdas et al.^l. We have translated 
the latter from the reported velocity data and into dimension- 
less units by using the known solvent properties, leading to 
£o = 0.095 pNs//im. In Fig. [2j also a result from simula- 
tions of monodisperse hard spheres by Carpen and BradyGSl 
is shown, highlighting that our choice of soft-sphere poten- 
tial and polydispersity are not crucial for the present discus- 
sion. While in both simulations, probe and host-liquid parti- 
cles were of the same size, 8 = a s /a = 1, the quoted values for 
the experiment in Ref. [12] are a s = 2.25 jum and a = 1.1 /um, 
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leading to a size ratio of 8 « 2.05. This difference explains a 
shift in force scales when comparing the two sets of data. 

Solid lines in Fig. [2] are fits using the schematic model. 
Fit parameters are given in Table [TJ additionally, under 
the constraint that the fits describe reasonably the known 
(simulation) or anticipated (experiment) linear response, we 
found (K||,k_l) = (0. 25, 0.0625) a /k B T for the simulation, and 
(K||,Kj_) = (0.046, 0.01 \5)a/k B T for the experimental data to 
give satisfying results. These are not the only parameter sets 
that give reasonable fits; on the schematic level, no more phys- 
ical significance can be attached to the choice of parameters. 

As is evident, the model reproduces the qualitative features 
of the available data, in particular the steep change of £ al- 
ready in the liquid around a threshold of O(20 k B T /a) for the 
simulation and O(200 k B T /a) for the experiment (larger due 
to the larger force required to pull free a larger probe). As 
pointed out earlier, 14 this threshold force is in the schematic 
model precisely defined as the force that is needed to pull the 
probe particle free even if the suspending host is glassy. Reas- 
suringly, microscopic calculations for the hard-sphere model 
yielded values comparable to the ones found in our simula- 
tion. From the schematic model, we get for the critical force 
just at the glass transition F^ x « 63 k B T/ a for the simulation, 
and F^ x w 800 k B T /a « 3pN for the experiment. At higher 
densities, even larger forces are needed to break the cages. 
According to MCT, the host system is in an ideal-glass state 
for densities (p > (p c « 0.595 in the simulation; hence, at van- 
ishing external force, the single-particle friction coefficient is 
infinite since the tracer remains localized in the glass. The 
simulation data for (p = 0.62 still show finite values of £ at 
all forces, which may partly be due to aging effects; it may 
also reflect deviations of the simulated glass from the ideal- 
ized MCT description. 

Before turning to a discussion of the large-force behavior, 
let us briefly discuss the steady-state probe velocities corre- 
sponding to Fig. [2] This is the quantity originally obtained in 
the experiment, and is shown in Fig. [3] The plot highlights the 
difficulty to stay within the linear-response regime in dense 
colloids. Velocities smaller than \0~ 3 k B T / a^o ~ 10~ 3 Dq/<2 
have to be measured, where Do ~ jum 2 / s is a typical short- time 
diffusion coefficient for a colloidal particle of radius a ~ 1 jum. 
The resulting velocities are on the order of a few particle di- 
ameters per week, which clearly is a challenge. Such small 
velocities pose the additional problem of ensuring a typical 
enough sample of the configuration space required in the en- 
semble average, a problem also for the simulations. 

In both Figs. [2] and [3] dashed lines represent the predictions 
of the schematic model when setting ju± = 0, i.e., ignoring the 
contribution of fluctuations perpendicular to the external force 
to the friction coefficient. Fits with somewhat different param- 
eters also setting v s 2 = have been presented earlier;^ they 
are qualitatively the same. The resulting curves still show the 




Fig. 3 Steady-state probe velocities (v)«> corresponding to Fig.[2| as 
a function of external force: computer- simulation (a) and 
experiment (b). Lines are the same schematic-model fits as in Fig. [2] 

generic features of the rapid decrease in friction in the vicin- 
ity of the force threshold. The most striking difference is that 
for ju± = 0, one obtains ^(F ex — » oo) = £ 0j the free-particle 
friction coefficient, as explained above. This behavior is not 
supported by the data, which clearly show a larger, density- 
dependent high-force plateau in ^(F^). 

At least within the schematic model, including fluctuations 
perpendicular to the direction of the external force restores 
the nontrivial high-force plateau. Here, § S j_(t) remains a real- 
valued positive function whose decay is fixed by the equi- 
librium short-time relaxation of the probe particle. Conse- 
quently, a finite contribution remains for even as F QX — » oo. 

A comparison between simulation and experiment regard- 
ing the high-force plateau in ^(F ex ) is instructive. Replotting 
the experimental data as in Fig. [2] reveals that the data is not 
yet in the limiting regime F QX — >> oo for all densities, but our 
schematic-model fits suggest a confident extrapolation. We 
obtain the values shown in Fig. [4] as symbols. The difference 
between the two data sets is the result of two physical differ- 
ences in the systems: the ratio 8 = a s /a and the influence of 
hydrodynamic interactions. To demonstrate this, we have in- 
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Fig. 4 High-force limit of the microscopic friction increment, 
A^oo/^o = C(^ ex ~^ °°)/Co _ 1> obtained from simulation (circles; 
size ratio 3 = a s /a= 1) and experiment (squares; 5 = 2.05), as a 
function of packing fraction (p. Lines are the low-density estimates 
(dotted) by Squires and Brady,^ corrected for packing effects by 
the equilibrium contact value from the BMCSL equation of state for 
these size ratios (solid lines). 




io" 2 io" 1 10° io 1 10" 1 10° io 1 

t [(ma/kT) m ] t 



Fig. 5 Probe-density correlation functions (|)| (t) from BD 
simulations at q = 3.85/a and (p = 0.55 (left) and from the 
corresponding schematic model fit (right). Top (bottom) panels 
show the real (imaginary) part of the correlators. For the simulation, 
aF GX /(k B T) = 1,5, 15,35, 100, and 250 (right to left). 



eluded in Fig. [4] as dotted lines the result obtained by Squires 
and Brady. 19 Evaluating the hydrodynamic equations for <p — >> 
0, they obtained A^o/^o = (l/4)q>(l +8) 2 . A first-order ac- 
count for packing effects in the host liquid is obtained by re- 
placing the packing fraction (p in this expression by cp&/(cp), 
where gd(<f>) is the probe-host-particle contact value. For the 
latter, we use the form corresponding to the BMCSL equation 
of state for hard- sphere mixtures due to Grundke and Hender- 
sonJ22in the limit of vanishing probe-particle density, gj((p) = 

( l_ (p + 3cp5/(l+5))/(l-(p) 2 + 2cp 2 (5/(l + 5)) 2 /(l-(p) 3 . 
Note that this expression differs slightly from the more com- 
mon Carnahan- Starling result for 8 = 1, used, e.g., by Carpen 
and Brady .ESI As shown in Fig. [4] this theory describes reason- 
ably well the experimental data when we set 8 = 2.05. The 
increase due to size-ratio effects alone is about a factor 4 at 

Table 1 Fit parameters for the MCT schematic model 



e Vs A* 

BDcp = 0.40 -0.98 9 0.45 

HS (p = 0.45 -0.87 9 0.59 

BDcp = 0.50 -0.54 9.75 0.62 

BDcp = 0.55 -0.23 11.25 0.74 

BDcp = 0.57 -0.185 14.25 0.78 

BD (p = 0.62 0.185 33 093 

expcp = 0.45 -0.17 4.5 14.3 

expcp = 0.50 -0.07 7.5 20 

expcp = 0.52 -0.03 15 25 

expcp = 0.53 -0.02 30 50 

exp <p = 0.55 -0.015 45 100 



(p ^ 0.5. The contact- value-corrected low-density expression 
for 8 = 1 exceeds the BD results; the difference might be due 
to the polydispersity and the slightly soft potential used in the 
simulations, as both factors reduce &/(cp). At (p < 0.2, the 
simulation data approaches the low-density asymptote. 

It is worth pointing out that the friction coefficient at small 
forces, F &x — >• 0, shows completely different scaling than the 
one for F ex — >> 00 ; while the latter is well described by hy- 
drodynamics corrected for pair-density effects, the former is 
governed by the approach of the equilibrium system to the 
glass transition, where ^(F ex — >> 0) — » 00 as (p — >> (p c from 
the liquid side. This is clearly demonstrated by the data 
shown in Fig. [2] We have checked that a t low densities, 
A£(F ex 0)/A£(F ex -> 00) = 2 as predictedESEni a i so holds 
for our simulations (up to about (p w 0.2). Clearly, at higher 
density A^(F ex — >> 0) rises much more dramatically than the 
F ex — >> 00 friction coefficient. 

The strength of the MCT model is that it can not only ex- 
plain the available data for the friction coefficient or steady- 
state velocity, but also the underlying probe dynamics in terms 
of the tagged-particle density correlation functions. We have 
obtained both the stationary and, for some values of F ex , also 
the transient correlation functions from the BD simulation; 
both are shown in the left-hand panels of Figs. [5] and [6] as 
continuous and dotted lines, respectively, for a wave num- 
ber qa = 3.85, packing fraction (p = 0.55 and various external 
forces. In Fig.[5j the alignment q \\ F ex has been chosen, while 
Fig. [6] shows q _L F QX . Even at the largest forces we investi- 
gate, the difference between stationary and transient correla- 
tion functions is not qualitative. In the direction of the force, 
the transient correlation function decays faster than the sta- 
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Fig. 6 Probe-density correlation functions ^j_(t) from BD 
simulations at q = 3.85/a and (p = 0.55 (left) and from the 
corresponding schematic model fit (right). Top (bottom) panels 
show real (imaginary) parts of the correlators for forces as in Fig. [5] 



tionary one. We interpret this as being due to a "wave front", 
i.e., a region of increased density building up in front of the 
pulled particle, which slows down the structural rearrange- 
ment in this direction once it is established in the steady state 
(see also the discussion in Ref. [19]). Just after applying the 
external force, this wave front has not been built up for short 
times, and this is probed by the transient correlation function. 

We now focus on the stationary correlation functions that 
are easier to obtain from the simulation, since the difference 
to the transient quantity is not important for our purposes. In 
mode-coupling theory based on integration through transients, 
the transient correlation functions are the central quantities, 
and the evolution towards the stationary functions needs to be 
understood separately.^ Thus, the schematic-model correla- 
tors that are shown in the right-hand panels should, strictly 
speaking, be interpreted as transient correlation functions. 

One recognizes a number of generic features in the correla- 
tion functions, that are present in the simulation and correctly 
reproduced by our schematic model. First, (|)| (t) is a com- 
plex quantity, and it shows pronounced oscillations at large 
F ex that are indicative of a probe moving with a finite velocity 
that is faster than the diffusive exploration of its configuration 
space. 14 For small external forces, the correlation functions 
stay close to the equilibrium tagged-particle correlator recov- 
ered for F QX — » 0, indicating the linear-response regime. As 
F QX is increased, the slow relaxation of the correlator is in- 
creasingly accelerated, which is the microscopic origin for the 
pronounced force thinning observed in Figs. [2] and [3] 

Taking q JL F ex , one instead recovers a real- valued correla- 
tion function. Figure [6] confirms this prediction for the BD 
computer simulation: the calculated imaginary part is zero 
within the noise level. Here, too, large forces lead to a short- 
ening of the equilibrium slow relaxation time, but no oscil- 
lations are seen, as the probe particle does not perform net 
motion in the direction perpendicular to the force. At very 
high F QX , the schematic model predicts a saturation, where 
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Fig. 7 Microscopic friction 1 + A^(F ex )/^o for the BD simulation at 
(p = 0.55. Circles: simulation data from Fig.[2j squares (triangles): 
contributions in Eq. ([T} coming only from the simulated q _L F ex 
(q || F ex ) correlation functions shown in Figs. |]and[6] Filled (open) 
symbols use the stationary (transient) correlation functions. 



approaches an exponential decay on the short-time dif- 
fusion scale. In the simulation, this is not (yet) clear. The 
integral over this remaining force-independent decay is, in the 
schematic model, the cause for the non- vanishing high-force 
friction increment over the solvent friction. The model thus 
predicts the large-force friction to be dominated by details of 
the short-time motion, but not by the slow structural relax- 
ation. The good quality of the low-density predictions for A^oo 
shown in Fig. [4] supports this. 

Further support for the schematic-model interpretation as- 
signing distinct roles to the principal directions q JL F ex and 
q || F ex at large external forces comes from the simulation data 
in Figs. [5] and [6] In the microscopic theory, Eq. ^ is replaced 
by a wave- vector dependent integral over § s ^(t)$ q (t), taking 
into account fluctuations in all directions relative to the force. 
In principle, this integral could be evaluated with simulation 
data for the (transient) correlation functions, without making 
use of the MCT equations of motion. It is a formidable task 
to sample enough correlation functions with sufficient statis- 
tics in the simulation. To give a qualitative impression, we 
only evaluate the integral over the simulated correlators for 
the two particular wave vectors shown above separately; this 
corresponds to the two contributions in the schematic Eq. ([7]). 
The host-liquid correlation function needed in this equation 
has also been evaluated at the same wave number, q = 3.85/a, 
in the BD simulation. The result is shown in Fig. J7] for an ex- 
emplary packing fraction. The integral using q \\ F ex (triangle 
symbols) gives ^ « at large forces, i.e. A^(F ex ) — ^ 0. Also 



in the BD simulation, the oscillations in ^(t) in the direction 
parallel to the external force cancel each other, and the main 
contribution to the friction at high external forces stems from 
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density fluctuations to wave vectors perpendicular to the force 
(squares in the figure). Our data for the transient correlation 
functions is not of sufficient quality to be used in the integra- 
tion, although this should strictly speaking be done. For the 
highest forces we simulated, we give a numerical estimate for 
the integral value according to Eq. ([7]) using the transient cor- 
relators (open symbols in the figure). The difference is not 
large (slightly negative values in some cases indicate the mag- 
nitude of the error in the integration procedure, due to cutting 
off at too early times). This further corroborates the point, that 
wave- vector contributions perpendicular to the external force 
are prominent in the large-force limit. 

5 Conclusions 

We have presented a schematic mode-coupling-theory anal- 
ysis of active nonlinear microrheology investigating the dy- 
namics of a probe particle pulled by a strong constant force 
through a glass-forming dense colloidal suspension. Ex- 
tending previous analysis, 14 our schematic model captures 
the qualitatively distinct behavior of probe-density correlation 
functions for fluctuations in the direction of and perpendicular 
to the external force. The distinct features - complex- valued 
and at high force oscillating correlation functions in the di- 
rection of the force, in contrast to force-thinning real-valued 
functions perpendicular to F QX - reflect the property that the 
probe-particle density distribution is biased in the direction of 
the force and attains a finite steady- state velocity super seeding 
the quiescent motion at high enough forces, while perpendic- 
ular to the force, no net motion occurs. These features are 
confirmed in Brownian dynamics simulations. 

Within the schematic model, contributions to the micro- 
scopic friction coefficient ^(.F ex ) coming from the different 
directions start being distinguishable upon entering the nonlin- 
ear regime. Both contributions drop sharply around a thresh- 
old that is large compared to forces induced by thermal fluctu- 
ations alone. They lead to pronounced force thinning and the 
eventual derealization of the driven probe even from a glassy 
surrounding matrix, at a critical force F c ex . The contribution to 
the friction from fluctuations parallel to the external force van- 
ishes as F ex — )> oo, but the contribution from the perpendicular 
directions remains finite, giving rise to a non-trivial high-force 
plateau above the Stokesian pure-solvent value. The existence 
of such a plateau is in qualitative agreement with both simu- 
lations and experiment, and also with predictions from low- 
density expansions. 

It remains difficult to compare the mechanism causing a 
nontrivial high-force plateau, A^oo > 0, in our present MCT 
approach with the one by Squires and Brady, 19 but the two 
descriptions share a number of apparent similarities. In Ref. 
IT9l it arises from a singular boundary layer, where advection 
and diffusion (otherwise negligible far from the probe at high 



force) compete close to the particle; transport of fluid par- 
ticles out of the boundary layer requires motion parallel to 
the probe surface, i.e., perpendicular to the force. In MCT, 
probe fluctuations transverse to the direction of the applied 
force need to be included. In either case, fast diffusive mo- 
tion of the host fluid particles, characterized by the short-time 
self-diffusivity 19 or 1^ in the schematic model, dominates the 
high-force friction. Structural correlations, which would be 
characterized by the much smaller long-time diffusivity re- 
spectively the slow relaxation time of the correlators in Figs. [5] 
and|6j dominate the quiescent friction. 
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